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As a result of the interplay between the intrinsic off-diagonal terms of the dipolar interaction 
and an applied transverse field ift, the diluted LiHoxYi_xF4 system at x> 0.5 is equivalent to a 
ferromagnet in a longitudinal random field (RF). At low Ht the quantum fluctuations between the 
Ising like doublet states are negligible, while the effective induced RF is appreciable. This results in a 
practically exact equivalence to the classical RF Ising model. By tuning ift, the applied longitudinal 
field, and the dilution, the Ising model can be realized in the presence of an effective RF, transverse 
field, and constant longitudinal field, all independently controlled. The experimental consequences 
for D — 1, 2, 3 dimensions are discussed. 



Since its seminal discussion by Imry and Mafl*], the 
Ising model in random longitudinal magnetic field was 
found to have many interesting realizations in nature, 
and has been a subject of intensive theoretical and ex- 
perimental investigation for both short-range and dipolar 
ferromagnetic (FM) interactions (see Refs. [H, 0,0, HI and 
references therein). In particular, experiments are mostly 
done on diluted antiferromagnets (DAFM) ^] , which were 
shown to be equivalent in their static critical behavior to 
the random field Ising model (RFIM)[i,0. In the DAFM 
the order parameter is given by the staggered magnetiza- 
tion, and a constant applied longitudinal field plays the 
role of a tunable random field (RF). However, the straight 
forward realization of the RFIM in a FM system has not 
found an experimental realization, due to the difficulty 
of applying a local random magnetic field. 

Interestingly, the diluted LiHoxYi_xF4 system pro- 
vides the possibility of producing effective local magnetic 
fields in the longitudinal direction, due to the interplay 
of the off-diagonal terms of the dipolar interactions and 
an applied constant magnetic field Ht in the transverse 
direction d, 0, Eli. Importantly, as is further ex- 
plained below, due to the strong hyperfine interactions 
in this system, for Ht < 0.5T the applied transverse field 
does generate an effective longitudinal RF due to fluctu- 
ations between each Ising ground state (GS) doublet and 
its corresponding excited state [TO, [HI, but the coupling 
between the two Ising doublet states is neglig ible[^. Fur- 
thermore, despite the presence of some correlations be- 
tween the random fields (RFs), and their origin in the 
long range dipolar interaction, we show below that on 
length scales larger than the inter-vacancy distance the 
behavior of the system is equivalent to that of the un- 
correlated RF model. Therefore, despite the presence 
of a finite Ht^ and finite correlations on short scales, 
LiHoxYi_xF4 at Ht < 0.5T is a perfect realization of the 
classical RFIM with a tunable random magnetic field. 

As Ht is enhanced, the effective transverse field be- 
comes appreciable. Thus, the Ising model Hamiltonian 
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can be realized in the LiHoxYi_xF4 with an effective lo- 
cal RF term 7^, an effective transverse (quantum) term 
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FIG. 1: Schematic picture of the phase diagram at some di- 
lution x< 1. At Ht < 0.5T the quantum fluctuations between 
the Ising doublet states are negligible, while the effective ran- 
dom field is already appreciable. In this regime the system 
is equivalent to the classical RFIM. At T < 0.2K only the 
electro- nuclear ground states are occupied, and the quantum 
phase transition in the presence of the random effective field 
can be studied. 



A, and an effective constant field in the longitudinal 
direction h\\. In the regime considered in this Rapid 
Communication, where x> 0.5, the system is FM at low 
T, Ht [H, [13] , and Jij correspond to the long range dipo- 
lar interaction. In comparison to the DAFM systems, 
the LiHoxYi_xF4 system is different since in the former 
the effective interaction is short range, and advantageous 
since in the DAFM system one can realize neither an ef- 
fective longitudinal field in the staggered magnetization, 
nor an appreciable quantum term. Importantly, the three 
effective fields. A, /ly and 7 = ^/Var{^ can be inde- 
pendently controlled by tuning Ht^ the applied magnetic 
field in the longitudinal direction Hz, and the dilution 
X. Specifically, for Ht <C IT and appreciable dilution 
7 > A, where A > 7 for i^t > 1-5T and x = 1 - x < 1. 

This gives for the first time an opportunity to study 
experimentally the RFIM in a FM system, in both the 
classical and quantum regimes. In particular, the realiza- 
tion of the Hamiltonian ([1]) in a FM system allows both 
the study of long-standing questions using direct bulk 
probes such as magnetization and susceptibility, and for 
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the first time an experimental study of models in which in 
addition to the random field a constant longitudinal field 
or a quantum term plays a role. Thus, the present work 
could be a basis for numerous new experimental possi- 
bilities, a few of which are the following: (i) The ther- 
mal and quantum phase transitions (PTs) of the RFIM 
at D = 3 can be studied in the whole T^Ht phase di- 
agram (see Figd]). Experiments studying the classical 
phase transition in 3D DAFM systems show that the RF 
prevents the sample from reaching equilibrium on exper- 
imentally accessible time scales |15l [iq . The study of 
the phase transition in diluted dipolar magnets can shed 
light on this interesting phenomenon, and the effect of 
the dipolar interaction on the theoretical predictions for 
the equivalence of the quantum and classical PTs in the 
RFIM [13, [3 can be tested as wen[i9|. Of particular 
interest is the ability, resulting from the novel quantum 
term, to control the relaxation rate of the system by tun- 
ing the effective transverse field, (ii) As a result of the 
possibility to apply tunable random field and constant 
field in the longitudinal direction, one could study ex- 
perimentally the hysteresis in the RFIM as function of 
Hz. Theoretical studies of this model in the absence [13 
and presence [21] of long range interactions have shown 
an intriguing disorder-driven out of equilibrium PT, a 
paradigm for the study of crackling noise in self-organized 
critical systems, (iii) In D = 1,2 dimensions, the predic- 
tions of Imry-Ma[l| and Binder j22| for the destruction of 
long-range FM order [2^ and the rounding of the phase 
transition in 2D can be verified by studying the magne- 
tization as function of the Ht dependent effective RF. 

All our analysis and results below can be easily gen- 
eralized to any dipolar Ising magnet, as is explained in 
detail in Ref.[lO|. In the following we concentrate on the 
LiIIOxYi_xF4 compound since this enables us to present 
quantitative results for this system which is of prime ex- 
perimental interest. 

Realization of the RFIM in the LiHoxYi_xF4. — 
Anisotropic dipolar magnets in general, and the 
LiHoxYi_xF4 compound in particular, are considered to 
be realizations of the Ising model. As a result of the 
crystal field anisotropy the GS is an Ising like doublet, 
and all but the longitudinal dipolar interactions are effec- 
tively reduced. The application of Ht induces quantum 
fluctuations (QF) between the Ising ground states (GSs), 
resulting in an effective transverse field Ising model 
(TFIM). For the undiluted system both the thermal and 
the quantum PT to the PM state were observed [2^. 
The dominant interaction in the LiHoxYi_xF4 system 

is dipolar[2l i^dip = E^,,c./3 ^^f 1^ ^^e pure 

LiHoF4, the off-diagonal terms of the dipolar interac- 
tion, e.g. JfJj^ are not only effectively reduced, but 
cancel by symmetry. However, at finite dilution, x< 1, 
this cancelation does not hold. Furthermore, in Ref. [TO] 
it was shown that in the presence of a constant Ht, the 
off-diagonal terms, despite being effectively reduced, are 
significant since their presence changes the symmetry of 
the system. The system is no longer symmetric under 
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FIG. 2: Y? jy? = (^/m/cVox)^ in three dimensions is plotted 
as function of = 4*L^ for dilutions x= 0.6 (blue triangles), 
0.96 (black squares), and 0.996 (red circles). L is the length of 
the cube in unit cells. Averaging is done over 10^ samples for 
each data point. Where not drawn, errorbars are smaller than 
the point size. Inset: \^\l\fN as function of x for N — 4000. 
|E| = tJo^VN for X <Cl, with r]o = 4.12. 

the transformation Sz —Sz, and an effective local lon- 
gitudinal magnetic RF hj is generated. In the regime 
Ht < l^o/(/iB^) 

i 

see Ref.flo!]. The random local energy is given by 7^ = 
jj^BhjS. Qq ~ lOK is the measure of the anisotropy, given 
by the excitation energy between each of the Ising GSs 
to its relevant excited state [lo|, and 5* ^ 5.5 is the mag- 
nitude of the effective Ho spin. 

Importantly, hj ([2|) is independent of the directions of 
the spins at sites i, and therefore equivalent for the spin- 
glass (SG) and FM regimes. However, the energy gain 
due to the effective RF does depend on the spin config- 
uration. Thus, while the randomness in the energy gain 
in the SG regime is a result of both the random posi- 
tion and the random orientation of the spins, in the FM 
regime only the former are random. Furthermore, the 
RFs hj are not uncorrelated. For x <C 1 the correlations 
of hj at distances smaller than the typical inter- vacancy 
distance, are large. However, as is shown explicitly be- 
low, for all dilutions, for distances larger than the inter- 
vacancy distance, the typical energy gain in flipping a 
domain Ejm = 7j obeys the Imry-Ma behavior. In 
particular Var{EiM) c<: where is the number of 
sites in a domain. 

Let us now calculate explicitly = (^.^. Vfj^ /Vq)'^ = 
{l/cVo)'^Ejj^, where c = 2/iBHtS'^ /Qq and Vq is the 
magnitude of the dipolar interaction at distance of one 
unit cell along the x direction. First, note that for a 
single vacancy z, as a consequence of lattice symmetry. 
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E,- = 0. This results in Y.ij V^f = Eij Vf, where 
i, j run over occupied sites and j over the unoccupied 
sites. Consider first the case of x <C 1. Assuming that 
the interactions causing the RF, V^Y ^ nearest neigh- 
bor only, then in leading order E'jj^ is proportional to the 
number of nearest neighbor vacancy pairs, and therefore 
to . In analogy to the case of dipolar interaction in a 
random system [26] it is expected that for the random in- 
teraction leading to the effective RF the dipolar and short 
range interactions would be equivalent. Thus, we expect 
Var{EiM) cx N for ah dilutions and Var{EiM) (x x'^N 
for X <C 1, where correlations of three or more vacant 
spins can be neglected. 

In Fig|2] we plot Var(S)/x^ for cubes in 3D as func- 
tion of N = for x= 0.6,0.96,0.996. Here and 
below lengths are given in terms of lattice unit cells. 
The dilutions were chosen to include the cases where 
X <C 1, and effective fields are strongly correlated on 
small distances, x= 0.6 where correlations due to the 
relative positions of three or more vacancies are signif- 
icant, and x= 0.96 in between. For each x we obtain 
a good linear fit, Var(S)/x^ = {r]{'x.))'^N^ which corre- 
sponds to Var{EiM) = (/?(x)cxVb)^A^. In the inset we 
plot ^ Var(S) /N as function of x for a 3D cube with 

= 4000, showing limx^o V = with r]o = 4.12. 

We therefore conclude that for distances larger than 
the inter- vacancy distance the system is equivalent to a 
RFIM, with an effective uncorrelated RF 
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Tjj being a random number with {rjj) = and Var{r]j) = 
7^^/x. Similar results (not shown) were obtained for D = 
1. Interestingly, in Eq.© is similar to the RF obtained 
in the SG regime [10], with x replacing x and the constant 
prefactor being different for the two regimes due to the 
different spin configuration. 

While our results are general to any dipolar Ising mag- 
net, we now use the specifics of the LiHoxYi_xF4 system 
to show that despite the presence of Ht the analogy to 
the classical RFIM is possible. The Ising like anisotropy 
is dictated by a crystal field Hamiltonian having large 
and terms[25]. In addition, a large (J^ + Jl) term 
results in a mixing of free ion states with AJ^ = ±4. 
The electronic GSs are an Ising like doublet, denoted 
|T),|i), belonging to the = (±7, ±3, Tl, T5) multi- 
plets. Although the full splitting of the 2J + 1 = 17 
states is 1] ~ 600K, the first excited state | F), belonging 
to the Jz = (6,2,-2,-6) multiplet is only = lOK 
higher than the GS doublet. The electronic Hamilto- 
nian allows a coupling between the electronic GSs which 
is second order in Ht. However, the low energy physics 
of the LiHoxYi_xF4 system is strongly affected by the 
strong hyperfine (hf) interaction between the Ho elec- 
tronic angular momentum and its nuclear spin \2^ . 
In particular, each electronic GS is split to 2/ + 1 = 8 
nearly equidistant states with Iz = —7/2, ...,7/2 and a 



separation of 205mK between them[28|. Thus, the rele- 
vant Ising like single spin GS doublet is of electronuclear 
type, having a definite and opposite nuclear spin states, 
i.e. a = 11,-7/2) and a = ||,7/2). The tunneling be- 
tween the states a and a requires the flipping of the nu- 
clear spin, resulting in a very weak coupling at low fields 
[less than lO-'^K at Ht < IT (Ref. [13)]. Only when | F) 
becomes well hybridized with the electronic GSs ||), ||), 
QF between the states a and a become appreciable. This 
occurs at Ht ~ Qo/{jiiBS) ~ 1.5T (see a detailed plot of 
A{Ht) in Fig 2 of Ref. 0). 

We first consider the low T low Ht (LTLH) regime 
{Ht < IT and T <C 0.2K). In this regime only the 
states a, a are appreciably occupied, and the quantum 
tunneling between these two states is negligible. Thus, 
the system is effectively a classical Ising model [9]. Im- 
portantly, the effective RF (|3]), originating from fluc- 
tuations between each Ising GS and its corresponding 
excited state[lo|, [HI, can already be made appreciable, 
since its behavior is linear in Ht. Thus, in this regime, 
which is of particular interest for the study of the de- 
struction of long range order (LRO) by the effective RF 
in = 1, 2 (see below), and the equilibrium and nonequi- 
librium phenomena in the ordered state in = 3, the 
analogy between the LiHoxYi_xF4 system and the RFIM 
is practically exact. 

Phase transitions in three dimensions. — As a result 
of the very different dependence of 7 and A on i^t, the 
thermal and quantum PTs can be well separated, as is 
sketched in FigHJ Both phase transitions occur outside 
the LTLH regime. For the classical PT one has to go 
to T > 0.2K, since x * 1.54K[T3|. Therefore, ah 

nuclear states are appreciably occupied near the transi- 
tion. Still, at low enough transverse field, Ht ^ 0.5T, 
QF are small between any of the time reversed states 
II, — T^), II, T^) [27]. Thus, the nuclear spin states can be 
traced over, and the classical PT is equivalent to that of 
the RFIM. Note, that the magnitude of the RF can be 
tuned by varying the dilution, and for each given sample 
by changing the magnitude of Ht, as can be inferred from 
Eq. (O. 

Considering the quantum phase transition (QPT) we 
focus on T <C 0.2K, where only the GS doublet is appre- 
ciably occupied. However, we still have to go out of the 
LTLH regime since strong QF are required. Specifically, 
in the dilutions of interest the QPT occurs at Ht > 3T. 
In Ref. [9] it was shown that QF between the electronu- 
clear Ising states become appreciable at Ht ~ 1.5T, and 
at Ht ^ 3T the GS has all hf levels well mixed [13 • How- 
ever, since at x> 0.5 the typical magnitude of the dipolar 
energy VqxS'^ ^ 0.2K, it is only when the QF of the elec- 
tronic spins are larger than the interaction that the PT 
to the PM phase takes place. At x= 1 it was argued that 
the system is equivalent to the TFIM, with the hf inter- 
actions effectively re- normalizing the transition field . 
For x< 1 and VqxS'^ ^ 0.2K we expect the same picture 
to be valid, with the addition of the effective longitudinal 
RF. 
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At the QPT Ht ^ Qo/il^BS), Eq. ([2]) is not valid, and 
QF between the Ising doublet states enhance the effective 
RF. The largest RFs that can be obtained are of the 
order of VoS/fiB- By changing the dilution x, one can 
therefore control the magnitude of the typical random 
magnetic field h at the QPT in the regime < h ^ 
Vb5'//iB, equivalent to < 7 ^ J, where J is the average 
nearest neighbor interaction. Thus, for x <C 1 the QPT 
can be studied in the presence of weak RF, while for a 
substantial dilution the RF, as well as the fiuctuations, 
are of the order of the interaction. Note, that LRO is 
destroyed when either QF or the RF are of the order of 
the interaction, and the situation where all three energy 
scales are of the same order is of both theoretical and 
experimental interest. 

Measuring the correlation length in one and two 
dimensions. — Two is the critical dimension for the sta- 
bility of the FM phase to the RFpj,i2j], as the energy gain 
and energy cost to flip a domain are both linear in the 
domain size. However, Binder [2§| has shown that due 
to the roughening of the domain wall, a two dimensional 
FM Ising system in a RF has a finite correlation length 
oc exp (Vb5'^/7)^. Here we consider first a strictly 2D 
system, where the plane includes the easy z-axis and is 
chosen as to have a FM phase at zero transverse field. We 
first consider an in-plane magnetic field transverse to 
the easy axis. For Hx <C flo/{jiiBS) the effective RF is 
given by Eq. (j3|), and therefore, for an appreciable dilu- 
tion, say 0.5 < X < 0.8, ex exp [((^q/mb^x)^]- Since 
5^1, > exp (5^) in the perturbative regime. As 
Hx is increased, QF between the Ising like doublet be- 
come significant, h increases and ^ decreases. Assuming 
one can reach the regime where <C the finite corre- 
lation length can be detected by studying M as function 
ofHx 



M ^ Mo{Hx)^^^ ^ Mo{Hx)^ , (4) 

as the ratio M/Mq depends linearly on ^. Here Mo{Hx) 
corresponds to the dependence of the magnetization on 
Hx in a FM domain due to single spin fiuctuations, as 
can be measured in a 3D sample where the system stays 
FM for the relevant magnetic fields [Hf <C flo/{jj.BS)]. 

Unlike systems such as the DAFM 
RB2CoxMgi_xF4[29], where the two dimensional- 
ity is due to the inherent layered structure, for 
the LiHoxYi_xF4 a realistic system would be quasi 
two-dimensional, with a finite width d. The energy 



cost to fiip a domain is then linear in d, while the 
energy gain is oc V^. As a result ^ is enhanced, 
^Q2D ^ exp [l]Q(i/(/iB^x)^], and therefore one has to 
keep d small in order to observe a finite ^. In comparison 
to layered systems, the observation of a finite ^ is more 
difficult. However, the study of the gradual cross-over 
between two and three dimensions is made possible. 
Note, that the degree of two dimensionality can be 
studied by measuring the dependence of the RF on 
the angle of Ht in the xy plane. This also allows the 
continuous change of the RF in a given sample. In 
particular, for the strict 2D case there is no RF when 
Ht is perpendicular to the plane. 

Using simple Imry-Ma[l] arguments one can show that 
for a quasi- ID system parallel to the z-axis and of cross- 
section A, the finite correlation length is given by ^ 
A{Qo/ jiiBHt)'^ ^ and the magnetization by 



«=.M„m)iy|=M„(fl,)/|-2|-. (5) 

Thus, the observation of M/Mq oc 1/Ht in the regime 
Qq^^/aJL <c; /iB^t <C ^o/S is the experimental mani- 
festation of the instability of the FM phase in one di- 
mensional Ising system to RF, as predicted by Imry and 
Ma. 

As explained in detail in Refs. [13, llH, our analysis 
above applies directly to a general anisotropic dipolar 
system, with a crystal field Hamiltonian e.g of the form 
HcF = DS^^ and with no hyperfine interactions. In- 
terestingly, in the LiHoxYi_xF4 system, if one neglects 
the hyperfine interactions, the peculiar form of the crys- 
tal field Hamiltonian which allows transitions between 
the two electronic ground states in second order per- 
turbation, leads to different physical results. Most im- 
portantly, the effective transverse term becomes appre- 
ciable at small i^t, and the RFIM is obtained at the 
expense of significant quantum fiuctuations. Thus, the 
proper consideration of the hyperfine interactions in the 
LiHoxYi_xF4 system is not only crucial for obtaining 
the correct effective fields, but also to obtain the quali- 
tative equivalence of the LiHoxYi_xF4 system to general 
diluted anisotropic dipolar magnets. 
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